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Abstract: It has been argued that the effect of cosmo logical structure formation 
on the average expansion rate is negligible, because the linear approximation to the 
metric remains applicable in the regime of non-linear density perturbations. We 
discuss why the arguments based on the linear theory are not valid. We emphasise 
the difference between Newtonian gravity and the weak field, small velocity limit of 
general relativity in the cosmological setting. 
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1. Introduction 

The backreaction conjecture and perturbation theory. It has been pro- 
posed that the observed increase in the expansion rate and the distance scale of the 
universe at late times relative to the matter-dominated homogeneous and isotropic 
Friedmann-Robertson- Walker (FRW) model could be explained by the breakdown 
of the homogeneous and isotropic approximation because of the formation of non- 
linear structures [1-5]. The effect of dumpiness on the average expansion rate is 
called backreaction [6-8]; see [9-11] for reviews. The exact Buchert equations for the 
average expansion rate show that large variance can lead to accelerated expansion as 
faster regions come to dominate the volume [8]. This effect has been demonstrated 
with exact toy models [10,12-16]. At late times there are deviations of order unity 
in the expansion rate between different regions, so this mechanism could also work in 
the real universe. The correct order of magnitude and timescale of the change of the 
expansion rate have been shown to emerge from the physics of structure formation in 
a semi-realistic model without any free parameters [17,18]. The relation between the 
average expansion rate and observations of light is also understood, though it should 
be established more rigorously and details remain to be worked out [19-21]. How- 
ever, there is no fully realistic calculation yet, and whether backreaction is important 
in the real universe remains an open question. The difference between Newtonian 
gravity and the weak field, small velocity limit of general relativity [22-29] plays an 
important part in the problem. Therefore, quantifying the importance of the growth 
of structures on the average expansion rate requires treating a statistically homoge- 
neous and isotropic but locally complicated non-linear system in general relativity. 
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However, it has been argued that the effect of non-hnear structures on the ex- 
pansion rate can be evaluated in hnear perturbation theory around the FRW met- 
ric [3,16,30-39], sidestepping subtleties of non-linear general relativity and the New- 
tonian limit. The argument is that even though the density perturbation becomes 
non-linear when structures form, the corresponding metric perturbation in the lon- 
gitudinal gauge, calculated from the Poisson equation, remains much smaller than 
unity, so the effect on the averages is negligible. There are multiple problems with 
this argument. Evaluating the effect on averages requires going at least to second 
order, so using first order perturbation theory is inconsistent, observables are not 
given by the metric alone, but by the metric and its derivatives (which can become 
large) and, finally, the linear equations do not, in fact, apply once the density per- 
turbation becomes non-linear. In short, it is not enough to calculate the magnitude 
of the effect in linear perturbation theory, the applicability of the linear treatment 
also has to be considered. 

Some of these arguments have been addressed before [5,8,10,17,21,40-42]. How- 
ever, as they are being repeated in the literature, it may be useful to discuss the issue 
in more detail than in [5,10,17,21], and from a slightly different perspective than 
in [8,40-42]. In section 2 wc consider perturbation theory around the FRW metric 
and show why the linear and second order calculations arc not sufficient for evaluat- 
ing backreaction once the density field becomes non-linear. Wc then look at the full 
non- linear equations for the averages and consider the Newtonian limit. In section 3 
we discuss previous work on this topic, and in section 4 we summarise our conclusions 
and outlook. 

2. Perturbations and the average expansion rate 
2.1 The perturbative calculation 

The Einstein equation and the metric. We assume that matter and geometry 
are related by the Einstein equation 



where is the Einstein tensor and Tq,^ is the energy-momentum tensor; we use 
units in which 87rGN = 1- where Gn is Newton's constant. We assume that the 
matter can be described as dust. 



where p is the energy density and u"' is the velocity of the observers, taken to be 
comoving with the dust. 



— , 



(2.1) 



Tap = pUaUp , 



(2.2) 
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The perturbed FRW metric in the longitudinal gauge is (we consider only a 
spatially flat background and only scalar perturbations) 

ds^ = -(1 + 2$(t, x))df + (1 - 2^{t, x)) aitfSijdx'dx^ . (2.3) 

The Einstein tensor for the metric (12.31) is 



G\ - -3H^{1 - 2$) - 2a-2V^* + QH^ 
G\ - -{2H + 3H^){1 - 2$) + 2^ + 6H^ + 2H^ 
-a-2v2(^ - $) + a-^dli^ - $) 

Go. ^ 2di{^ + H^) , (2.4) 

where ~ denotes dropping terms which are higher than first order (or, later, second 
order; it should be clear from the context which is meant) in $ or \E', dot denotes 
derivative with respect to the background coordinate time t, and no summation is 
implied in We also split the velocity into the background and the perturbation, 
= + and assume that 6u°' is small. From the normalisation condition 
gapu°'u^ = — 1 it then follows that m° ~ 1 — $. 



With the energy-momentum tensor ( p.2|) and the Einstein tensor (2^), the Ein- 
stein equation ( p.l|) reduces, at first order, to 



2,H^{1 - 2$) + 2a-2v2* - 6i/^ ~ p (2.5) 
2H + 3i7^ - 2^ - QHi! - 2H^ + a"V^(^ - $) - a'^dli^ - $) = (2.6) 

didji^ -<!?)= (2.7) 
2di{4f + H^) ~ -p6ui . (2.8) 

Note that we have not made any assumptions about the perturbations of p. From 
(13) it follows that ^ - $ = A{t,x^) + B{t,x^) + C{t,x^), where A, B and C are 
arbitrary functions. We are mostly interested in the situation when the perturbations 
are statistically homogeneous and isotropic, in which case A = B = C = 0, and we 
assume this from now on. (The condition ^ — $ = would also follow from the 
technical requirement that the Fourier transform of \E' — $ exists.) 

The static case. Let us first consider the static case H = 0, and choose a = 1. 
The set of equations (^3| ) -( ^^) reduces to 

2V^$ ~ p (2.9) 
$ = (2.10) 
2di^ - -p5ui . (2.11) 

To be consistent with neglecting terms which are second order in $, we should 
discard the right-hand side of ( |2.11| ), because according to ( p^.9| ), p is of order $. We 
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then obtain the result $ = At + B[x), where A is a constant and B[x) is determined 
by the density via ( p.9|) . It is possible for p to have large variations without $ 
becoming large or the first order treatment becoming invalid. (This is the case in 
the solar system, for example). However, in that case there is a slight inconsistency 
in the treatment, because we have assumed that = |p is small. If p is allowed 
to be large, we should equally treat as a large term, so products such as <I>V^$ 
should not be discarded. However, we should then take into account second order 
terms in the metric, because they can be of the same order. Let us look at this in 
more detail in the cosmological situation. 

The cosmological case. With H ^ 0, the Einstein equation ( p.5|) -( p^ reads 

3i/^(l - 2$) + 2a~2V^$ - 6//$ ~ p (2.12) 
2ij + 3/f2 - 2$ - 8//$ ~ (2.13) 
2di{^ + //$) ~ ~p6ui . (2.14) 

As is usual, we assume that the background and first order equations are sepa- 
rately satisfied. This follows if we assume that the average of $ over the background 
space vanishes. We split the density into the background value and the perturbation, 
p = p + 6p, but do not assume that 6p is small. We then have 

3H^ = p (2.15) 

2a-^V^^ -QH^^ -QH^ = 6p (2.16) 

2H + 3H^ = (2.17) 

(i> + 4:H^ = (2.18) 

Sui = --di{^ + H<^) . (2.19) 
P 



Equations ( p. 17] ) and (|2.18|) , which come from the pressure-free condition (p.l3|) . 



determine the evolution of a and $, regardless of the energy density. They lead to 
the standard relations a oc t^/^, $ = A(x) + B(x)t^^^^, where A and B are arbitrary 
functions. According to (|2.15|) and (|2.16|) , the density contrast 6 = Sp/p is related 
to $ by 

' = 3(^^'* - 2* - I* . (2.20) 

We have nowhere required that 6 should be small, so one could at first sight 
think that ( p.20| ) applies even when 6 becomes of order unity, as long as $ remains 
small, analogously to the static case. (In the static case p = 0, so 5 is not defined, 
but the variation of p between different regions of space can be large.) Keeping to the 
linear theory, this is not true. The time evolution of $ is determined independently 
of p, and inserting $ = A + Bt^^/^ into (^.2(j|) shows that, dropping the decaying 
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mode, the density contrast has a constant part and a part which is proportional to a. 
Therefore S grows without hmit. For an underdense region, S cannot go below —1, so 
this evolution is clearly not correct as 6 becomes of order unity. The behaviour is also 
wrong for overdense regions, as is well known from the spherical collapse model [43] 
(see [44] for reviews). 

For the volume expansion rate 6 = Vq,m" we have 



where we have on the second line dropped the decaying mode of $. (For the ex- 
pression in terms of the proper time measured by the observers, see [5].) It is clear 
that the expansion rate given by the linear theory is wrong when 6 becomes of order 
unity. 

The reason for the breakdown. We have assumed that the Einstein tensor and 
the velocity can be expanded linearly in the metric perturbations. We have found 
that the observables calculated using this procedure fail to describe the real behaviour 
when 6 ~ V^$/(aif)^ becomes of order ±1, even if $ would seem to remain small, 
so the linearly perturbed metric would appear to be valid. 

The reason is that in neglecting all terms which are second order in $, we have 
implicitly assumed that terms such as $V^$/(aif)^ are much smaller than i.e. 
that |V^$/(aif)^| ~ \6\ ^ 1. To extend the calculation into the regime \6\ > 1, 
we would have to expand to second order in $. But to do this consistently, we 
have to include the intrinsic second order terms in addition to the squares of first 
order terms. Indeed, the distinction between the two is gauge- dependent [45] , as first 
order quantities are not invariant under second order gauge transformations. And at 
second order, the metric cannot be written in the simple diagonal form (|2.3| ) [46] . It 
may be that the effect of the second order terms is small for a particular quantity of 
interest, but this has to be determined via a consistent calculation. 

Effectively, there are three expansion parameters: $, di^/ (aH) and V^$/(ai7)^ ~ 
S. (No higher derivatives appear, because the Einstein equation is second order.) 
The formal perturbation expansion is defined in powers of the metric perturbation, 
treated as an infinitesimal quantity [47]. However, in the real universe, the metric 
perturbation has a finite amplitude, so the gradients can make the other expansion 
parameters large even when $ remains small. The gradient is a dimensional quantity, 
so a comparison scale must enter. In cosmology, the relevant scale is aH, and since 
aH decreases in a decelerating FRW universe, gradients become more important with 
time. We can also view this as follows: for a time-independent $, the magnitude 
of is fixed in time, while the curvature scale of the universe, to which it is 
compared, decreases. This kind of an instability is not present in the static case. 



e ~ 3i7 - 3(<l> + H<^) + diu' 
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In a situation with multiple expansion parameters, perturbation theory can be 
expected to remain valid when all parameters are small, and to fail when all of 
them are large. When some parameters become large while others remain small, 
the validity of perturbation theory depends on the system, and on the quantity 
under consideration. In cosmology, the metric ( p.3| ) is simply the first order term in 
an expansion, and when gradients of the metric perturbation become large, higher 
order terms can no longer be neglected. This does not necessarily mean that all first 
order results are wrong: a consistent calculation with higher order terms may show 
that some linear relations are valid. However, this cannot be determined using the 
linear theory. 

For backreaction, the important quantity is the average expansion rate. (The 
primary quantities are of course observables defined in terms of measurements of 
light; for the connection to the average expansion rate, see [19-21].) One might 
argue that even if the linear theory fails to correctly describe the local quantities 
when density perturbations are non-linear, the effect on the averages nevertheless 
remains small. To address this issue, let us see what happens when we expand to 
second order. 



The average expansion rate at second order. Taking the metric ( |2.^j| ) and 

calculating 9 = Vq,u° to second order in $, we obtain (we take $ = 0; see [5] for the 
general expression) 

where Hr = 2/ (3r) is the background expansion rate in terms of the proper time r 
of comoving observers. The last term, with four derivatives, is of order 5^, so there 
are large local variations in the expansion rate. Averaging ( |2.22 ) on the hypersurface 
of constant proper time, we obtain [5] 

/ 22 1 22 1 

^(9.(V2<|.9.$))o) , (2.23) 



27 (aH) 

where () is a proper average with the correct volume element, ()o is an average taken 
on the background hypersurface of constant proper time, without perturbations in 
the volume element, and we have assumed ($)o = 0. It is noteworthy that the term 
with four derivatives, which has the largest amplitude locally, is a boundary term. 
Before discussing this feature, let us note that this calculation is not consistent, 
because we have used the first order metric to calculate a second order quantity, i.e. 
we have neglected intrinsic second order terms. To obtain a result which does not 
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depend on the gauge, it is necessary to truncate the metric consistently at second 
order instead of first order. The result is then [45]^ 



dinal gauge happens to give qualitatively the right answer, but the coefficients of the 
terms are wrong. (In first order perturbation theory, doing the calculation in the syn- 
chronous comoving gauge, for example, would give a qualitatively different result.) 
Note that there is nothing in the result of the first order calculation that would indi- 
cate that the answer is wrong. An average of a total derivative can be converted into 
a surface integral of a flux through the boundary. If the distribution is statistically 
homogeneous and isotropic, there is no preferred direction, so the integral vanishes 
(up to statistical fluctuations). (In perturbation theory, the technical requirement 
that $ can be expanded in Fourier modes would lead to the same conclusion.) With 
vanishing boundary terms, the correction to the mean is ~ {di^di^)o/{aH)^, which 
is of the order 10~^ for a realistic linear theory power spectrum. To see the failure of 
the linear theory expanded to second order, we have to work with the second order 
metric. With the metric truncated at first order, it is impossible to determine the 
magnitude of the higher order terms which are neglected. As the intrinsic second 
order terms are as large as the first order terms squared, the question arises as to 
the magnitude of the terms which are even higher order. While calculating the coef- 
ficients of the various terms would be an involved task, it is straightforward to write 
down their general form. 

The general structure of the corrections. At second order, the possible correc- 
tion terms are the squares of the three expansion parameters, ($^)o, {di^di^)o/ {aHY 
and (V2$V2<l>)o/(ai^)^ ~ (5^)0. (The quantity ($V2$)o is equal to -{d.^di^)^ up 
to a boundary term.) For simplicity, we take $ = 0. As long as |$| ~ taking 
into account time-dependence would simply introduce more terms of the same order 
of magnitude, and would not lead to any qualitative change. When the average ex- 
pansion rate is expressed in terms of the proper time r (as opposed to the unphysical 
coordinate time t), $ appears in the expansion rate only with derivatives acting on 
it [48]. This is to be expected, because if $ depends only on time, it corresponds to 

-•^As an aside, in [45] it is assumed, as is usual, that the equations are satisfied separately order 
by order in perturbation theory. While this procedure is self-consistent, there seems to be no 
rigorous justification for it beyond first order. At first order, the equations for the background 
and perturbations decouple, assuming that the average of the perturbations vanishes. Starting at 
second order, the average of the perturbations does not vanish, so decoupling of the background 
and perturbations is an extra assumption. 





Comparing ( p.23| ) and ( p.24| ) shows that the first order calculation in the longitu- 
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using a different time coordinate, not to having a physical degree of freedom. Assum- 
ing that the higher order equations are satisfied order by order and the perturbations 
are Gaussian, all higher order terms from scalar perturbations factorise into products 
of these three expectation values, since they are sourced by the first order terms. In 
contrast, vector and tensor perturbations (which necessarily arise at higher orders) 
have solutions which do not need to be supported by a source, so their contribution 
cannot be completely expressed in terms of the first order seed fields. (See [46] for 
the second order case.) However, such terms are expected to be subdominant to the 
scalar perturbations in the non-linear regime. 

We now return to the feature that at second order, the term with the highest 
number of derivatives (and therefore locally the largest amplitude) is a boundary 
term and as such vanishes upon averaging, up to statistical fluctuations. In [5] it 
was argued that at higher orders there might not be such a cancellation for the 
leading terms. However, in [49] it was realised that because each factor of di/ (aH) is 
accompanied by one power of the speed of light c, the terms with the highest number 
of spatial derivatives are the ones which dominate in the Newtonian limit c — ?■ cxo. 
In Newtonian gravity, the backreaction correction is exactly a boundary term [7]. 
Thus, in general relativity the term with the highest number of spatial derivatives 
at each order in perturbation theory is a boundary term^. The general structure of 
the corrections from scalar perturbations to the average expansion rate is therefore 
(see also [40]) 

(^)=3i/. (^l + ^(9.$9.<l>)ogA„(52)^ + ...j , (2.25) 

where A„ are constants and . . . indicates subleading terms with a smaller number 
of derivatives, such as ((9j$9j$)™(5^)Q~"'/(aif)^'", with n > m > 1. In powers of 
$, the Xn term is of order 2n + 2. The term (5^)o can appear at fourth order in 
$ at the earliest, where the leading correction is (9j$(?i$)o/(aif )^ x {6^)o. This 
term grows without bound with increasing \6\, so the breakdown of the perturbative 
expansion is transparent. None of the coefficients A„ have been calculated. It is 
possible to determine Ai in third order perturbation theory, which is being developed 
[50], without a full fourth order calculation [51]. However, calculating Ai would be 
inconclusive, because at every order, there are an increasing number of terms that 
grow even faster as |5| becomes of order unity. 

In [40] it was argued that the series ( |2.25D would have only a finite number of 
gradient terms when $ is taken as the full metric perturbation and not only the 
linear part (and ^ 7^ $ is included). However, this is not the case: for a metric of 



^Assuming that the leading order general relativity result reduces to the Newtonian theory at 
all orders. As we discuss in section 2.2, this is not necessarily true. If that is not the case, the series 
(2.25) is even more divergent. 
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the form (p.3|) , the velocity (and thus also 6) expanded as a series in $ necessarily 
contains an infinite (or zero) number of spatial derivatives [5]. And as we have noted, 
beyond first order, the metric cannot be written in the form ( p.3| ) [46]. 

From the fact that the series ( p.25| ) would naively seem to diverge at (5^)o ~ 1 we 
cannot conclude that the sum of the correction terms would be large. However, we 
can definitely say that the series expansion does not prove that the correction would 
be small when $ is small. The magnitude of the effect has to be established with non- 
perturbative methods, or a resummation of the series. For studies in the spherically 
symmetric situation where the exact solution is known, see [16,52-55]. In particular, 
[55] shows that it is possible to have a large effect on the observable distance-redshift 
relation even when the metric can be written in the form (|2.3| ) (at least on the 
lightcone). These models are not conclusive of the cosmological situation, which is 
not spherically symmetric^. Different resummation schemes have been applied in 
Newtonian cosmology [56], and it would be interesting if such methods could be 
extended to general relativity. 

We would still expect to recover linear equations for perturbations with wave- 
lengths much larger than the size of the structures, as is usual in statistical physics. 
These should look similar to perturbation equations around the FRW universe, with 
correction terms due to the underlying structure [10]. This is also suggested by the 
success of FRW perturbation theory in describing observations of large-scale struc- 
ture. Such equations would be analogous to the Buchert equations, which look like 
FRW equations with correction terms, though their physical content is different, as 
they involve only average quantities and not local expansion. The effect of backreac- 
tion cannot be described merely as a change in the FRW background [10,17,19,21], 
unlike argued in [38,57,58]. Even though the average expansion rate will always agree 
with that of some FRW model, other observables will in general not be the same as 
in that FRW universe. In particular, the relationship between the average expansion 
rate and the luminosity distance is different than in FRW models if backreaction is 
important [19,21,59,60]. 

We have discussed perturbation theory as it is most commonly formulated, by 
adding perturbations on top of a background (and previous perturbations). The al- 
ternative is to take the full non-linear system and linearise it. In cosmology (unlike 
in the spherically symmetric case) we cannot write down the exact solution to lin- 
earise. However, it is at least possible to build perturbation theory by starting from 
the full exact equations, written in the covariant formalism, and linearise around 
the FRW solution [61-63]. This has the benefit that all terms are included to begin 
with, so it is transparent to estimate what is being dropped, unlike in the case when 
perturbations are added to a background order by order. In addition, the covari- 
ant formalism deals only with measurable quantities and the physical spacetime, so 

^Note that in [54] the average spatial curvature is small, so it is clear that backreaction is not 
important. 
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there are no gauge artifacts. Instead of a perturbative analysis, we go directly to 
the exact non-perturbative equations for physical insight into the effect of perturba- 
tions becoming large. A comparison of the general relativistic and Newtonian cases 
is also instructive, given that backreaction vanishes in the latter, for a statistically 
homogeneous and isotropic distribution. 



2.2 The Newtonian limit 

The Buchert equations. If the matter is irrotational dust, the exact equations 
which describe the effect of inhomogeneities on the average expansion rate {9) = 
Sdrtt/a in general relativity are [8] (for the case with non-dust matter or rotation, 
see [21,64-66]) 

3^ = -47rGN(p) + Q (2.26) 
a 

3^-^ = 8nGM-le^R)-lQ (2.27) 

dr{p) + 3 — {p) = , (2.28) 
a 

where ^^^R is the spatial curvature and Q is the backreaction variable defined as 

Q^lm-{ef)-2{a') , (2.29) 



where cr^ is the shear scalar. The integrability condition between (|2.26|) and ( p. 27] ) 



is 



d^{(^)R) + 2^((3)/?) = -d,Q - 6^Q . (2.30) 
a a 

If Q = 0, we have {^'^^R) oc a ^; in particular, this holds for all exactly homogeneous 
and isotropic universes [67]. The system of equations ( p.26|) -( p.28D closes once we 
are given Q or {^^^R); because of the integrability condition ( ^.30[ ), the effect of 
dumpiness can be viewed equivalently in terms of either quantity. The Raychaudhuri 



equation ( ^^.26] ) together with ( ^^.29] ) shows that, apart from a possible {^^^R) oc a" 



2 



term, deviations from homogeneity and isotropy have a large effect on the average 
expansion rate only when the variance of the expansion rate is large, and is not 
cancelled by the shear (or the shear is large, and is not cancelled by the variance). 
This shows that (5^)o = 1 is not a sufficient condition for a large effect on the average 
expansion rate. However, it is necessary that the deviation of the expansion rate from 
the mean is large in a large fraction of space (assuming that the deviation is at most 
of the same order of magnitude as the mean, which is true in cosmology). 
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Newtonian gravity and beyond. In Newtonian gravity, tlie counterparts of tlie 
Raycliaudliuri equation (|2.26|) and tlie conservation equation ( p.28|) are identical to 



the relativistic equations, but there is no analogue of the Hamiltonian constraint 
( p.27|) . The variance of the expansion rate and the shear combine to give a total 
derivative, so Q reduces to a boundary term [7]. (If the vorticity is non-zero, it is 
included in this boundary term.) Thus, if backreaction is important in a statistically 
homogeneous and isotropic universe, this must be due to non-Newtonian aspects of 
general relativity [10,11,17,40,41,49,68]. 

In the expansion (|2.25|) , all correction terms are post-Newtonian. The term 



c^{di^di^)o/{aHy, which may be identified as the square of a pecuhar velocity, 
v'^/c^, suppresses the post-Newtonian terms. However, the terms it multiplies can 
become very large as c^(V^$V^$)o/(a-ff)^ ~ ('^^)o grows. This demonstrates that 
in general relativity, non-Newtonian effects can be important even when velocities 
are small and fields are weak. An exact example is given by rotating and expanding 
dust. In general relativity, there are no dust solutions which have non-zero expansion 
and rotation but zero shear [69]. However, in Newtonian gravity such solutions are 
exactly known [22,25]. Analysis of the Newtonian theory in this case would be 
misleading, because the Newtonian solutions betray no sign of the fact that starting 
from general relativity, they do not exist, even at small velocities and weak fields^. 

This issue arises due to the indeterminacy of Newtonian cosmology, which is 
related to the absence of Newtonian analogues of the magnetic component of the Weyl 
tensor^ and the evolution equation of the electric component of the Weyl tensor [22, 
24,26,28,29,71,72]. Newtonian cosmology is only defined up to boundary conditions 
[24,28,29]. This shows up in the fact that Q is a boundary term, and the average 
expansion rate is determined by what happens at the boundary. In general relativity 
this is not the case, and backreaction is given by integrals over the volume. 

If the volume considered has periodic boundary conditions or is statistically 
homogeneous and isotropic, then in Newtonian gravity Q vanishes, and the first 
integral of the Raychaudhuri equation ( p. 261) leads to an equation which looks like the 



Hamiltonian constraint ( p.27| ) with Q = and {^^^R) = Ea~'^, where E is a constant 
of motion which may be identified with (being proportional to) the conserved energy 
of the isolated system. In relativistic cosmology, the conserved energy is replaced 
by spatial curvature, which has no physical analogue in Newtonian gravity, so the 
interpretation of this term is different in the two theories, even in the FRW case. 
There is no conservation law for the spatial curvature in general relativity, so {^^^R) 
can evolve in a non-trivial manner, unlike the total energy of an isolated Newtonian 
system [17,68]. The equivalent statement in terms of the backreaction variable Q 



■*This underlines that it is not sufficient to look at the Newtonian limit of the equations of general 
relativity, but it is necessary to consider the limit of solutions, because in general the operations of 
taking the limit and solving the equations do not commute [23,25,70]. 

^We can equivalently say that in Newtonian gravity the magnetic component vanishes identically. 
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is that in Newtonian gravity the variance of the expansion rate always equals three 
times average shear scalar (up to a boundary term), as we see from ( p.29|) , while in 
general relativity there is no such constraint. 

In second order perturbation theory we have, dropping boundary terms and 
taking 6 = 0, 

c;n 1 

^ {dm^)o , (2.31) 
9 

where the second approximation holds if $ is small compared to its gradient. Second 
order relativistic perturbation theory around a spatially flat FRW background is 
close to Newtonian gravity in the sense that there is an exact cancellation between 
the variance and the shear, so Q = 0. The variance and the shear can be calculated 
using first order theory, because they vanish for the background [51]. There is only 
a single non-Newtonian term, the spatial curvature, which is proportional to 
according to the integrability condition (|2.30|) . To determine the coefficient of this 



term, it is necessary to go to second order. 

Already in first and second order perturbation theory the variance and the shear 
are large. A large backreaction effect in general relativity does not require the vari- 
ance of the expansion rate to be larger than expected, it is enough that the cancel- 
lation with the shear is not perfect, unlike in Newtonian cosmology. However, the 
spatial curvature does have to become large [10,17,73]. In a realistic cosmological 
setting, this is easy to understand. The spatial curvature of the initial overdense 
and underdense regions averages to zero in the linear regime, but once perturbations 
become non-linear, the evolution of overdense and underdense regions is different, 
and the average will in general deviate from zero. It is to be expected that if the 
volume of the universe becomes dominated by underdense voids which expand faster 
than overdense regions, the average spatial curvature will be negative. 

Comparing Newtonian gravity and general relativity in cosmology is different 
than in the case of isolated, asymptotically fiat systems. For isolated systems, both 
Newtonian gravity and general relativity are well-defined. In contrast, while rela- 
tivistic cosmology is well-defined there is no unique Newtonian theory of cosmology, 
because the Newtonian equations are only defined up to boundary terms which have 
to be specified at all times [24,28,29]. (We could alternatively say that Newtonian 
gravity is a theory of isolated systems only, and there are an infinite number of 
possible generalisations to the cosmological setting.) This shortcoming of the New- 
tonian theory is often hidden in cosmology by the assumption of periodic boundary 
conditions (sometimes implicitly through the use of Fourier series). For periodic 
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boundary conditions, Newtonian cosmology does have a well-defined initial value 
problem [28,74], at least to all orders in perturbation theory, but this situation does 
not correspond to the real universe, which is not periodic on the observed scales. 

So in the cosmological context, post-Newtonian corrections refer to the difference 
between a realistic relativistic cosmological model, which is thus far not tractable, 
and a Newtonian model which is defined only when it does not correspond to the 
real situation. In particular, it is not possible to estimate backreaction from usual N- 
body simulations, because the relativistic degrees of freedom are absent. Relativistic 
cosmological simulations would in principle provide an answer. The initial condi- 
tions, the matter model and the equations of motion are known and well-defined, 
but solving the system in full generality is not computationally feasible. It would 
be interesting to obtain a reduced system that would retain the relevant relativistic 
degrees of freedom while being tractable. Because relativistic cosmology is not sensi- 
tive to boundary conditions in the same way as Newtonian cosmology, the periodicity 
required for a numerical implementation would not be a crucial limiting factor. 

One check on the correctness of the Newtonian treatment in the non- linear regime 
is provided by comparison of N-body simulations with observations of structures. (In 
the simulations, matter with negative pressure or modified gravity is introduced to 
change the background expansion rate. Without such an addition, observations of 
the expansion rate and the distance scale are, of course, already completely dis- 
crepant with the Newtonian model.) On large scales significant differences between 
the simulations and observations have been reported. The observed homogeneity 
scale is an order of magnitude larger than in N-body simulations [75], and the num- 
ber of very luminous superclusters is about five times larger in observations than in 
simulations [76]. Whether this reflects a deficiency of the Newtonian treatment, or 
instead indicates a problem in the way simulations arc done or the observational data 
is analysed is not clear. The discrepancy could also be due to an incorrect choice of 
initial conditions, matter content or theory of gravity. 

Relativistic dust models which are Newtonian-like are a very restricted class [77] . 
The close relation of linearly perturbed relativistic FRW models and Newtonian 
gravity may be misleading because Newtonian-like models suffer from a linearisation 
instability. In general relativity, the Newtonian constraint that the magnetic com- 
ponent of the Weyl tensor vanishes identically is, in general, not propagated in time. 
However, in the linear theory, the constraint is trivially satisfied at all times. There 
are thus linear theory dust solutions which are not the limit of any non-linear solu- 
tion. More importantly, this shows that relativistic dust models do not, in general, 
have Newtonian counterparts and their evolution cannot be described in Newtonian 
theory. For practical applications in cosmology, the important issue is the quantita- 
tive importance of the non- Newtonian features, which depends on the solution under 
consideration. For addressing this question it would be useful to understand better 
the relation between the evolution of the electric and magnetic components of the 
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Weyl tensor and spatial curvature in the context of cosmological structure formation. 



3. Comparison to previous work 

Arguments in linear theory. There have been various claims that backreaction 
in the real universe is negligible [3,16,30-39]. In particular, it has been argued that 
the relative magnitude of backreaction corrections is given by the square of the pe- 
culiar velocity. All these studies, except for [31,38], expand quantities calculated 
with the first order metric to second order, which, as we have seen, is not in general 
consistent. Often the physical expansion rate, proper time and hypersurface of av- 
eraging are also not correctly identified. That was first done in [5], while the first 
consistent second order calculation was done in [45]. Let us discuss some of the other 
shortcomings of these studies. 

A numerical estimate of the correction to the expansion rate from expanding the 
first order metric to second order was first given in [30]. In [31] the Zel'dovich approx- 
imation was used to obtain the second order metric and calculate corrections to the 
expansion rate. In [3] the correction to the Einstein equation was calculated in the 
same manner. In [34,36,38] a calculation similar to the one in [5,45] was done, with 
some variations. (In [36] , the perturbation of the volume element was inconsistently 
neglected.) In [32] a similar calculation for the correction to the 00-component of the 
Einstein equation was done, and non-linear scaling relations were used for the density 
power spectrum, but this cannot compensate for using only the linear metric. As dis- 
cussed in section pTl], such calculations lead to the correction term (5j$9j$)o/(a-ff)^, 



which is of the order 10^^. (In [31], the value 10~^ was obtained instead.) This result 
is the origin of the idea that the magnitude of backreaction is given by the square of 
the peculiar velocity, because at second order we have {di^di^)o/ {aHy ~ {uiU'')o. 
However, beyond second order, the expansion parameter ((5^)o is also involved, so the 
second order calculation is inconclusive, and ((9i$9j$)o/(a-ff)^ -C 1 is not a sufficient 
condition for small backreaction. 

As an aside, note that is the (spatial component of the) deviation of the phys- 
ical velocity of observers comoving with the dust fluid from a fictitious background 
velocity. This is a coordinate-dependent quantity, and we can always set = by 
choosing coordinates which are comoving with the observers. In order to determine 
a physical peculiar velocity, we have to define another physical velocity field to com- 
pare u"' to [78]. (In the longitudinal gauge in the linear theory, UiU^ does give the 
physical magnitude of the deviation from uniform motion.) 

In [33] it was asserted that the linear metric ( p. 3D (or the equivalent with a 
spatially curved background) describes the universe on all scales, except in the vicin- 
ity of black holes and neutron stars. It was then claimed that if the conditions 
|<I>| < 1, < a^'^di^di^, (di^di^Y < didj^didj^ hold, non-linear corrections 
are negligible. In fact, the metric (|2.3[) cannot (with a dust source) simultaneously 
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describe the static metric in the solar system and cosmological expansion, as is clear 
from the expressions for 9 and 5 in (|2.20| ) and ( p.21|) together with ( p.l8|) . We can 
either use the metric with if = to describe a static structure or take if 7^ to de- 
scribe a cosmologically evolving region, but these conditions are obviously mutually 
incompatible. Apart from the inconsistency of using the linear theory to calculate 
second order quantities, if we nevertheless took ( p^.3D and expanded observables in 
terms of then at fourth order and higher we would expect to obtain terms involving 
(5^)0, which are not necessarily small. (Of course, these corrections are meaningless 
without accounting for the intrinsic higher order terms.) 

It was also argued in [33] that since the average expansion rate depends on the 
choice of the averaging hypersurface, accelerated expansion could arise as a gauge 
artifact. However, we should distinguish three different concepts, namely gauge de- 
pendence, coordinate dependence and dependence on the averaging hypersurface. 
Gauge dependence arises due to ambiguity in the mapping between the perturbed 
physical spacetime and a fictitious background spacetime. When points of the fic- 
titious and real spacetime with the same coordinate values are taken to map to 
each other, gauge dependence reduces to choice of coordinates, but in general it is 
a distinct issue. In the covariant formalism with the full non-linear equations, we 
deal only with physical quantities and the real spacetime, so there is no gauge issue. 
As all quantities are defined covariantly, independent of the choice of coordinates, 
the dependence on the coordinate system appears only in the usual transformation 
properties of tensors under coordinate changes. In particular, covariantly defined 
averages of scalar quantities such as the volume expansion rate do not depend on 
the coordinate system when expressed in terms of a physical observable such as the 
observer's proper time [45]. They do, however, depend on the choice on the averaging 
hypersurface [48,79,80]. The reason is that the averaging hypersurface is physical 
issue, unlike coordinates or gauge. For irrotational dust, there is a preferred folia- 
tion which is orthogonal to the fluid flow, and which coincides with the hypersurface 
of constant proper time [10,17,41]. However, the hypersurface should be chosen 
based on analysis of observables, and cannot be determined on abstract mathemati- 
cal grounds [10, 19,21]. Any average quantities are of course useful only insofar they 
give an approximate description of what is actually measured. (For discussion of 
gauge-invariance in averaging, see also [66,81].) 

In [35] the linear metric was again expanded to second order. It was assumed that 
the average energy density is the same as the background energy density, which is not 
true beyond first order. Accordingly, one obtains equations which are inconsistent 
[17]. It was also argued that backreaction vanishes in a 2-1-1- dimensional model. This 
is not surprising, because in 2-1-1 dimensions, the integral of the spatial curvature is a 
topological invariant. Therefore it cannot evolve in time, similarly to the total energy 
of a Newtonian universe discussed in section |27^ [17]. Therefore, in 2-1-1 dimensions, 
backreaction can only give Q oc a~^, and it is not clear whether the coefficient can 
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be non-zero. (This follows from the 2+1- dimensional analogue of (|2.30|) : in d + 1 
dimensions, the last term on the right-hand side is —2d^Q.) 

In [16,37] an iterative calculation was done in the macroscopic gravity formalism 
[82], which is an extension of general relativity. We are interested in what happens 
in general relativity, but let us note that the method of [16,37] was to take the first 
order metric, expand to second order to obtain a new background and then repeat 
the process. This way, one never moves beyond first order perturbation theory. As 
we have seen, in general relativity the higher order terms in general do not have the 
same form as the first and second order terms, so this kind of an analysis would not 
be correct there. 

In [39] quantities were expanded to second order in the linear metric, with the 
usual result. The correction (9i$(9j$)o/(aiJ)^ xp ~ v'^p was identified with a pressure 
term. For clarity we note that the pressure measured by observers comoving with 
dust is zero by definition. The physical interpretation of the second order correction 
is spatial curvature, not pressure. To determine the pressure (and anisotropic stress 
and energy flux) generated in the process of structure formation, it is necessary to 
go beyond the ideal fluid treatment [83]. For the importance of non-dust terms for 
backreaction, see [21]. 

Let us also comment on some studies which claim not a small, but instead a 
possibly large backreaction effect from perturbation theory. 

A series expansion similar to (|2.25| ) was presented in [49]. It was argued that 
the expansion parameter is not (5^)o, but a quantity which becomes of order unity 
around the present time. However, the expansion in [49] is incorrect, because it does 
not take into account the factorisation of higher order terms into two-point functions 
and the constraint that a non-zero two-point function must contain an even number 
of momenta [17]. The only preferred era in the perturbative expansion is (5^)o = 1, 
signifying the formation of the first generation of gravitationally bound objects. For 
typical models of supersymmetric dark matter this happens around a redshift of 40- 
60 [84] , considerably earlier than the present day. (As discussed above, the failure of 
the simple perturbative expansion ( p.25|) at (5^)o = 1 does not alone indicate that 



backreaction would be large.) 

In [85] first order theory expanded to second order was used to estimate backre- 
action and compared to observations. Apart from the question of the applicability of 
first (or second) order perturbation theory, the correction term used is qualitatively 
wrong, because the momentum scale in the integral is misidentified with the size of 
the averaging domain (the suppression of the leading correction due to the fact that 
it is a boundary term is also neglected); see section 5.1 of [17]. 

4. Conclusion 

Summary and discussion. If we consider cosmological perturbations which are 
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initially small and Gaussian with zero mean, it is necessary to go at least to second 
order to find their effect on the average expansion of the universe, called backreaction. 
Nevertheless, taking the first order metric in the longitudinal gauge and expanding 
quantities to second order serendipitously gives almost the correct second order result 
[5,45]. At second order, the perturbations only lead to a spatial curvature term 
with an amplitude of 10~^. Several papers have thus claimed that backreaction is 
negligible in the real universe, based on the argument that the linear metric is a good 
approximation even when density perturbations are non-linear. 

We have discussed why the linearly perturbed FRW metric does not in general 
correctly describe the situation once density perturbations become non-linear. The 
effect on the average expansion rate vanishes at the linear level by construction, and 
at second order, the intrinsic second order terms are of the same order as squares 
of the first order terms (in fact, the division between the two is gauge- dependent). 
At higher orders, generic correction terms become larger than unity as density per- 
turbations become non-linear. This does not necessarily mean that the effect on the 
average expansion rate is large, simply that the naive perturbative expansion is no 
longer valid. 

The important question is not in which form the metric can be written, but 
what happens to measurable quantities. For this purpose it is useful to consider the 
covariant formalism, which deals only with physical degrees of freedom and is fully 
non-linear. The effect of deviations from homogeneity and isotropy is quantified by 
the Buchert equations, which show that the average expansion rate will significantly 
deviate from the FRW behaviour when the variance of the expansion rate is of order 
unity and does not cancel against the shear (or vice versa) [8]. Even in the lin- 
ear (and second order) theory, the variance of the expansion rate becomes of order 
unity as density perturbations become non-linear. However, there is no significant 
backreaction in relativistic second order theory because the variance cancels exactly 
against the shear apart from a boundary term, a feature shared by non-linear Newto- 
nian gravity. In exact general relativity, there is no such cancellation, so Newtonian 
theory is not sufficient for evaluating backreaction. 

Determining whether structure formation in the real universe leads to a large 
enough variance for backreaction to be important requires dealing with a locally 
complex non-perturbative system in general relativity. However, details of the lo- 
cal behaviour are not needed, statistical information about the distribution of the 
expansion rate in different regions is enough. A semi-realistic statistical calculation 
found a rise of 10-30% in the expansion rate relative to the FRW value around a 
time of 10 billion years [17, 18], which agrees with the observations within an order 
of magnitude. The calculation involved several approximations, and a more careful 
treatment is needed. In particular, the difference between relativistic and Newtonian 
cosmology should be better understood to isolate the relevant relativistic degrees of 
freedom, related to spatial curvature and the electric and magnetic components of 
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the Weyl tensor. 
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